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Q^' 1. Introduction 

Oh. The idea of parton saturation naturally arises in deep inelastic scatter- 

(— I I ing (DIS) at small values of the Bjorken variable x, studied currently at the 

HERA collider in DESY. The main measured small-x effect is a strong rise of 
the proton structure function F2 in the limit x —>■ 0, for fixed virtuality Q^. 
^ ■ Interpreted with the help of linear evolution equations of Quantum Chromo- 

dynamics, this is a reflection of increasing parton densities (sea quarks and 
gluons). Formally, the parton rise is so strong that for sufficiently small x 
the computed cross section violates unitarity, which signals that important 
physical effects were neglected in the approximation leading to the linear 
evolution equations. These effects are responsible for saturation of the par- 
ton densities by taming their strong rise. Whether the saturation effects 
are already visible in the F2 at HERA is disputable, the linear DGLAP 
evolution equations describe the data for > 2 GeV^. However, more ex- 
clusive diffractive processes strongly hint towards the affirmative answer to 
the posed question. Moreover, the transition region of ~ I GeV^ for F2, 
where the DGLAP analysis encounters basic difficulties, is easily described 
within an approach based on parton saturation. An important result of 
saturation is the existence of a saturation scale which is reflected in a new 
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scaling law for inclusive DIS cross section. The small-x data do really show 
such a regularity. 

In the following we present basic concepts leading to the notion of parton 
saturation and discuss the relevance of this effect for experimental data, 
mainly from HERA. 



2. Linear evolution equations 

In the standard description of ep DIS, the structure function F2 is de- 
termined by the quark and antiquark densities in the proton, 

F2(x,g2) = ^e^x{g/(x,Q2) + g/(x,Q2)} + o{as), (1) 
/ 

where the summation is done over quark flavours. As usual, x = Q^/(2P-g) 
and = —q^, where P and q are the proton and virtual photon momenta. 
The quark densities together with the gluon density g{x,Q'^) satisfies the 
DGLAP evolution equations [1]. In the matrix notation, q = {qf,qf,g), 



/ -P{x/z,Q')ci{z,Q^) (2) 



dlnQ^ 

where the matrix of splitting functions is computed perturbatively, 

P(z,Q2) = as{Q^) P«(z) + a2(g2) p(2)(^) ^ 

Eq. (2) allows to determine the scale dependence in (1) provided initial 
conditions at some scale Qq are given. This is done by fitting their form in 
X to the existing DIS data. In this way logarithmic scaling violation of F2 
is explained. The above description applies to ^ ^qcd perturbative 
QCD to be valid. How low in and in x the pure DGLAP approach is 
applicable should be determined from the analysis of data. In principle, the 
values of x and indicating that boundary are correlated, i.e. = Q^{x). 

The low region brings an issue of the (l/Q^)" corrections (twist ex- 
pansion) to F2, treated systematically in the operator product expansion 
of two electromagnetic currents. Eq. (1) takes into account only the loga- 
rithmic contribution in this expansion. For moderate x this is a dominant 
contribution, as shown by successful DGLAP analyses of DIS data. For 
X <^ 1, however, each term in the twist expansion, 

F2{x,Q^) = (x,lnQ2) + Fi'\x,lnQ^)^ + (4) 

is equally important due to the presence of large logarithms ln(l/x). Thus, 
a new systematics is necessary allowing to resum these large logarithms in 
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Fig. 1. k_\_- factorization in DIS, eq. (5). 



the small- a; limit, independent of the twist expansion. This is done through 
the A;j_-factorization formula [2] which contains all twists 



where $(Q^,k) is the virtual photon impact factor describing the process 
7* — > — >■ 7*, see Fig. 1. The function /(x,k) is the unintegrated gluon 
distribution, related to the gluon density at large and small x by 



In the above k is a two-dimensional vector of transverse momenta of two 
gluons which couple to the quark loop in four possible ways. Notice the 
integration over all values of |k| in (5). This means that the region of 
small momenta needs special attention to avoid problems with perturbative 
stability. The unintegrated gluon distribution satisfies the BFKL equation 
[3] which can be written as an evolution equation in the rapidity Y = 
ln(l/x). In the leading order approximation, when the terms proportional 
to (q!s ln(l/x))" are resummed, 



Diagrammatically, the BFKL equation resums the ladder diagrams with 
two exchanged (rcggcizcd) gluons interacting through the real gluons in the 
rungs strongly ordered in rapidity. Each emitted gluon leads to the large 
logarithm a5ln(l/x). Since this is a colourless exchange which gives the 
dominant behaviour at large energy -y/i of the 7*p system {x = Q'^/s), it is 
termed the BFKL pomeron. 
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Y=ln(l/x) 




Fig. 2. Diffusion pattern from the BFKL equation. 



The solution of eq. (7) gives the x dependence of F2 at small x. In the 
limit X ^ and fixed, it is given by the azimuthally symmetric saddle 
point solution 



where D ^ ln(l/x). Thus for decreasing x, the solution features a power- 
like rise with x and diffusion in Infe^. Even though at some initial x = xq 
the solution is concentrated around k'^ ^ kg, it diffuses into the region of 
larger and smaller values of k"^ for a; ^ 0, see Fig. 2. Diffusion into the 
infrared region means that the BFKL approximation is strongly sensitive 
to nonperturbative effects since the cross section might be dominated by 
a contribution from |k| ~ A. This problem can be phenomenologically 
cured by a separate treatment of the infrared domain [4] . However, infrared 
diffusion signals fundamental difficulty which reflects incompleteness of the 
BFKL approximation. This is intimately related to violation of unitarity 
reflected in the Froissart bound: F2 < chi^(l/.x), in contrast to the power- 
like rise in (8). The next-to- leading corrections to the BFKL equation, see 
[5] and references therein, weaken the rise but the two discussed problems 
remain. 

The missing physical effect in the BFKL approximation is parton satura- 
tion. In the infinite momentum frame of the proton, the following picture of 
this phenomenon exists. The transverse size of gluons with transverse mo- 
mentum k is proportional to l/|k|. For large |k|, the BFKL mechanism of 
gluon radiation, g gg, populates the transverse space with large number 
(per unit of rapidity) of small size gluons. The same mechanism also applies 
to large size gluons with small transverse momenta. In this case, however. 



f{x,k^ 
Vk^ 



) 



= X 



12a. In2/. exp(- ln^(fcVfcg)/Z^) 



(8) 
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Fig. 3. The illustration of factorization in eq. (10). 



the BFKL approach is incomplete since large gluons strongly overlap and 
fusion processes, gg g, are equally important. By taking these processes 
into account unitarity is restored and infrared diffusion cured. The key ele- 
ment for this is the emergence of a rising with rapidity (energy) saturation 
scale Qs{Y), which is a fundamental property of parton saturation [6, 7]. 
The density of large gluons with |k| < QsiY) no longer strongly rises and 
unitarity is restored. The low momentum contribution to the cross section is 
dominated by |k| ~ Qs{Y) 3> A.qcd for sufficiently large Y, which allows to 
avoid the problem of strong sensitivity to the infrared domain. In the forth- 
coming we show how this idea is realized in practice. For complementary 
introductions, see [8]. 

3. The dipole picture 

The problem of restoration of unitarity in the description of DIS at 
small X was pioneered by Gribov, Levin and Ryskin in [6]. In there non- 
linear corrections to the DGLAP evolution equations were presented and 
subsequently rigorously derived in the double logarithmic approximation 
by Mueller and Qiu in [9]. In this approximation, is large and x is 
small such that terms proportional to (a^ InxlnQ^)" ~ 1 are resummed. 
We want to avoid the high assumption and approach the problem in 
the approximation in which the BFKL formulation of DIS was discussed. 
To this end we transform the A;^-factorization formula (5) into the Fourier 
conjugate representation where the transverse momentum k is traded for 
its conjugate transverse separation r. Thus the following formula is found, 
see e.g. [10], 

F2 = QV{4.7r^aem)i<JT + (TL), (9) 

and 

fTT,L(x,Q2) = J (fr dz |^'T,L(g',r,z)|2 a{r,x) (10) 
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Fig. 4. Emission of a gluon in the large limit. 

where (Jt,l are 7*p cross sections for the indicated polarizations. Here r 
is the qq dipole transverse separation and z is the longitudinal momentum 
fraction of the dipole quark/antiquark with respect to the momentum q' = 
q+xp. The physical interpretation of factorization in (10) is provided in the 
proton rest frame. The virtual photon 7* splits into a qq pair long before 
the interaction with the proton since the formation time, Tqq ^ l/(xM), is 
very large in the small-x limit. This process is described by the lowest Fock 
state light-cone wave function of the photon, ^t.l [H]- The interaction 
with the proton is encoded in the qq dipole cross section a{r,x). 

The BFKL approximation has a nice formulation in the dipole picture in 
the large Nc limit [12, 13]. In this limit, emission of a gluon by a quark can 
be viewed as a creation of two dipoles out of the parent dipole, see Fig. 4. 
Each new dipole can radiate a gluon with much softer longitudinal momenta, 
which leads to new dipoles. In this way the parent qq dipole r evolves into 
a collection of dipoles r' with the density n{r,r' ,x). The interaction of each 
dipole with the proton occurs through exchange of a single perturbative 
gluon. Thus in this picture the dipole cross section equals 

a{r,x) = j (fr'n{r,r',x)agp{r') (11) 

where agp{r') describes the interaction of a dipole r' with the proton The 
BFKL effects are located in evolution of the photon wave function, the 
large logarithms ln(l/a;) come from strong ordering of longitudinal momenta 
of the emitted gluons. Consequently, the dipole density obeys the BFKL 
equation in the r-space 

= !^[fjLUr + ,y,.) - e(r^-p^)„(r,r'.x)}. (12) 

oY vr J Trp^ J 

The variable r' is a parameter in eq. (12), thus the dipole cross section (11) 
obeys the same equation. Notice that the integral kernel on the r.h.s. of 
eq. (12) has mathematically the same form as in the BFKL equation in 
the momentum space (7). Like in the momentum space, the saddle point 



dis3 printed on February 1, 2008 



7 



solution exhibits the power- hke rise in x and diffusion in Inr^: 

n(r,r^.T)) _ ^_4iVea.in2/. exp(-ln^(rVr^2)/£>) 
vr \'kL) 

and D ~ ln(l/x). Now, infrared diffusion corresponds to an unhmited rise 
of the dipole density for large size parent dipoles. From (11) and (13), the 
BFKL dipole cross section behaves as for any r which leads to violation 
of unitarity when a; — as in the momentum space. Now, the saturation 
effects would correspond to interactions between dipoles. Before discussing 
QCD formulations of these effects, we present a phenomenological model 
which captures essential features of parton saturation and describes the 
DIS data. 



4. The saturation model 

In this model [14], the dipole cross section is bounded by an energy 
independent value ctq which assures unitarity of F2, 



(T(r,a;) = o"o s 1 — exp 



where T{jq{x)^ called saturation radius, is given by 

Rq{x) = 1 GeV"^ (x/xo)^/^ 



(14) 



(15) 



The three parameters, ao = 23 mb, A ~ 0.3 and xq = 3 - 10~^, were fitted to 
all small-a; DIS data with x < 10~^ . The resulting curves for different values 
of X are shown in Fig. 5. At small r, a features colour transparency, o" ~ r^. 
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Fig. 5. The dipole cross section in the saturation model. 



8 



DIS3 PRINTED ON FEBRUARY 1, 2008 



which is perturbative QCD phenomenon, while for large r, a saturates, a ~ 
(jQ. The transition between the two regimes is governed by the increasing 
with X saturation radius Rj^ix). Therefore, an important feature of the 
model is that for decreasing x, the dipole cross section saturates for smaller 
dipole sizes, see Fig. 5. The BFKL rise with x is encoded in the small r part 
of a since for r <C Rq, o" ~ x~^. In contrast to the BFKL result, however, 
with decreasing x the saturation radius Ro{x) gets smaller and consequently 
a never exceeds ao- This is the way infrared diffusion, which concerns large 
dipoles, is tamed by the existence of the saturation scale: 

Qlix)^l/Rlix) ^ x-\ (16) 

Notice that the saturation scale rises when x ^ 0. With the dipole cross 
section (14) several remarkable features of DIS data are described. 




Fig. 6. The cross section a^»p = ctt + for fixed 7*p center- of-mass energy W 
(x = Q'^/W^). The dotted lines, for which niq = 0, show the effect of the dipole 
quark mass niq = 140 MeV. 
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Transition to low 

First, the transition to low values in the proton structure function. 
This has an intuitive physical interpretation if we relate the saturation ra- 
dius to the mean transverse distance between partons. If the partonic sys- 
tem is dilute and the qq dipole probe with a characteristic size 1/Q \s much 
smaller than the saturation radius, 1/Q <^ Ro{x)j the logarithmic behaviour 
is found [14, 15] 

F2 - j^^ln{Q'Rlix)) . (17) 

In the opposite case, when 1/Q :$> Ro{x), the parton system looks dense for 
the probe and 

^^-^'^''-Iq^w)- '''' 

This transition is shown in Fig. 6 for a^^p = ar + (Tl ~ F2/Q^- With a 
nonzero quark mass in the qq dipole even the photoproduction data can be 
described, which reflects a kind of continuity in the proposed description. 
The transition region between the two regimes, defined by the condition 
(critical line in the (x, Q^)-plane) 

Q'Bl{x) = l, (19) 

is found to be around Q^ = 1 GeV^ at HERA kinematics (the solid line 
across the model curves in Fig. 6). The critical line also indicates the limit 
of vaHdity of the twist expansion (4) in the {x , Q^)-plaiie [16]. It is also 
interesting to analyze the energy dependence resulting from the saturation 
model. If this is done for increasing Q^ values, smooth transition between 
the soft {F2 ~ x^^'^^) and hard pomeron (F2 ~ x~^'^) values is found [14]. 



4.2. DIS diffraction 

Secondly, once the dipole cross section is determined from the inclusive 
data analysis, it can be used as a prediction for diffractive DIS [17] since for 
the qq pair in the diffractive final state 



dt I t=o 



= / d^'^dz |5'T,L(r,z)|2 a\x,r). (20) 



The total diffractive cross section a!^^p = a^l+aj^ is found after dividing (20) 
by the diffractive slope Bd, taken from the experiment. A very important 
result of the saturation model is that with the form (14), the constant ratio 
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as a function of and x of the diffractive and inclusive cross sections is 
naturally explained [15], 



a. 



D 

7*p 



1 



£7- 



7*p 



HQ^Rlix)) 



(21) 



With the diffractive qq and qqg components, the description of the diffractive 
data is quite satisfactory, see Fig. 7 [17]. 

The diffractive interactions are very important for tracing saturation 
effects since they are mainly sensitive to intermediate dipole sizes, r > 2/Q, 
which directly probe the saturation part of the dipole cross section. In the 
inclusive diffractive cross section, in contrast to the fully inclusive case, the 
r < 2/Q contribution is suppressed by additional power of [17, 15]. 
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Fig. 7. The diffractive data from ZEUS [18] against the saturation model curves. 
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4-3. Geometric scaling 
Thirdly, the particular feature of (14), i.e. its scaling form, 

a{r,x) = a{r/Ro{x)) , (22) 

leads to the observation [19] that the virtual photon-proton cross section 
(or F2/Q'^) is a function of only one dimensionless variable r = Q^Rq{x) 
instead of x and separately 

(^^*p{x, Q"^) = o-^*p(r) . (23) 

The new scaling (called geometric scaling), valid in the small x domain, 
is shown in Fig. 8 for the DIS data [19]. The scaling variable r is the 
ratio of two scales, the photon virtuality and the saturation scale (16). 
Therefore, in its essence geometric scaling is a manifestation of the existence 
of the saturation scale. 




Fig. 8. Geometric scaling for small-x DIS data with a; < 10 ^, the scaling variable 
r = Q^Rl{x). 
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It is interesting to look at the Regge limit, x —>■ and fixed. From 
(18), the saturation model leads to 

F2 ~ Q^ln{l/x) (24) 

which is in accord with the Froissart bound, straightforwardly applied to 
DIS: F2 < cln^(l/x). The logarithm in (24) comes from the wave function 
ill eq. (10), thus there is still room for additional logarithmic dependence 
on X of the dipole cross section. 



5. Nonlinear evolution equations 

The phenomenological success of the model (14) raises the question 
about its relation to QCD. In particular, to what extent is the transi- 
tion to saturation justified by perturbative QCD? The saturation scale 
Q'^{x) ~ 1 GeV^ for x 10^^, thus we hope that at least the onset of 
saturation could be described by weakly coupled QCD. 

The problem of saturation effects was attacked by Kovchegov in the 
dipole picture [20]. He found in the large Nc limit a nonlinear equation 
which is a special case of more general hierarchy of equations derived by 
Balitsky [21]. Strictly speaking, the derived equation is only justified for a 
scattering on a large nucleus since a certain class of diagrams, which could 
be relevant for the proton, is suppressed for DIS on a nucleus. However, this 
equation can be considered as a QCD model for ep DIS at small x. The basic 
quantity in Kovchegov's formulation is the forward scattering amplitude of 
the qq dipole (originated from the virtual photon) on a nucleus, A/"(r,b,F). 
If X and y are the transverse positions of the dipole quarks with respect to 
the center of the nucleus then r = x — y and b = (x + y)/2, where the latter 
quantity is the impact parameter of the dipole. The dipole cross section is 
given by 

a{r,x) = 2 j (fhN{r,h,Y). (25) 

Comparing eq. (25) with eq. (14), we see that in the saturation model 
the following hypothesis on the form of N is made 



iV(r,b,y) = 1-exp -— _ 9(6 -60) (26) 




where ctq = 27r6g. The theta function can also be shifted into the argument 
of the exponent. Physically, the saturation model corresponds to the proton 
being a disk in the transverse plane with a sharp boundary. Any edge 
effects are ignored. Saturation leads to a uniform blackening of the disk 
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Fig. 9. Saturation (blackening) in the saturation model. Black means a = Uq. 

with decreasing x (for the parent dipole of a fixed size r) without changing 
the disk size, see Fig. 9. 

In the approach of Kovchcgov, the ampUtudc iV(r, b, Y) is determined 
dynamicahy from the nonhnear evolution equation derived in the large 
limit. Transforming back into the momentum representation, 



where = Ncag/n and K is the kernel of the linear BFKL equation [20]. In 
the dipole picture in the leading ln(l/x) approximation, the above equation 
describes not only the production but also merging of dipoles in the dipole 
cascade. In addition, an energy independent initial condition for the evolu- 
tion is supplied at Y = Yq, which describes the scattering of the parent qq 
dipole off the nuclcons through multiple two-gluon colour singlet exchange 
[22]. The effect of the nonlinear evolution equation is the generation of the 
interacting dipole cascade which brings the energy dependence of the solu- 
tion. In terms of Feynman diagrams, eq. (28) resums fan diagrams with the 
BFKL ladders and triple pomeron couplings in the large Nc limit [23], see 
Fig. 10. 

Assuming the sharp disk model for the 6-dependence, both analytical 
and numerical studies of eq. (28) [20, 23, 24, 25, 26] confirm the picture of 
blackening, anticipated in the saturation model. The dipole cross section 
(25) computed after solving eq. (28) features colour transparency for small 
r and saturation at large r with a similar dependence on x, i.e. for smaller x 
the saturation occurs for smaller dipole sizes. Moreover, the saturation scale 
emerges in such analyses. This is illustrated in detail in Fig. 11 where the 
solution ^(fc, Y) of eq. (28) (solid lines) and the solution of the corresponding 
linear BFKL equation (dotted lines), projected on the (A;, y)-plane, are 




(27) 




(28) 



14 



DIS3 PRINTED ON FEBRUARY 1, 2008 




Fig. 10. The BFKL fan diagrams resummed by eq. (28). 

shown [26]. The BFKL solution exhibits diffusion into small and large 
values of k even though the initial condition is concentrated sA k = k^. The 
nonlinear saturation effects in the equation (28) tame infrared diffusion, 
pulling the solution into the infrared safe region. The straight lines for 
the nonlinear solution reflect the scaling property: (j){k^Y) = <j){k/Qs{Y)), 
where the saturation scale QsiY) corresponds to the central line k = Qs{Y), 
see [26] for more details. It is important to note that the scaling sets in 
independent of initial conditions and is strictly valid for k < Qs{Y). For 
k > Qs(Y), however, the scaling is mildly violated. This was analyzed 
recently in [27] where the region of validity of geometric scaling was found. 

Unfortunately, the agreement between the results based on the equation 
(28) and the saturation model is only qualitative and detailed studies of the 
DIS data with this equation are necessary. 

An intriguing question, related to proton confinement, is the impact 
parameter dependence of the dipole-proton forward scattering amplitude 
A^(r, b, Y). Certainly, the sharp disk picture with a fixed boundary is over- 
simplified if not questionable. Recently, detailed studies with eq. (28) in the 
r-representation were done [28], assuming the initial condition A^(r, b, Yq) = 
N{r) S(h) with some nonperturbative profile function S(h) which falls as 
exp{—2m-,rb) at large b. The result of this analysis is shown in Fig. 12. 
The black disk corresponds to saturation but now the saturation area rises 
according to the nonlinear evolution equation (28) (or the equation of the 
Colour Glass Condensate [29]) when x ^ 0. Thus for each impact parame- 
ter b unitarity limit is achieved. The grey area corresponds to lower parton 
densities where the linear BFKL equation applies. Within this approach 
the Proissart bound is saturated: F2 ~ ln^(l/x). 

The impact parameter dependence of the forward scattering amplitude 
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logi„(k/lGeV) 

Fig. 11. Diffusion in the BFKL and Balitsky-Kovchegov equations. 



can be studied through the t-dependence of difFractive vector meson pro- 
duction in DIS [30], see Fig. 13. In this model, the vector meson production 
is a three step process, the virtual photon splits into the qq pair which in- 
teracts with the proton with the amplitude ^gf ~^, and then forms a meson 
described by the wave function 




Fig. 12. Black disk expansion and saturation. 



16 



DIS3 PRINTED ON FEBRUARY 1, 2008 




Fig. 13. Diffractive vector meson production. 

where My is the vector meson mass. The amphtude 

Al'l-P{r,A,x) = 2 J d%N{r,h,x)e^-^ (30) 

and A is a two dimensional vector of transverse momentum transferred from 
the proton into the qq system, t = — A^. The presented above results are 
confirmed qualitatively although the proton is not fully black at the central 
impact parameter 6 = in HERA kinematics. We would like to stress that 
diffractive vector meson production in scattering is potentially the best 
process to study saturation effects in DIS since the transverse size of the qq 
pair forming a meson is controlled by the vector meson mass 

r = 1/^M2 +g2. (31) 

Thus we expect saturation effects to be more important for larger (lighter) 
vector mesons. The first analysis in this direction, done in [31] with the 
saturation model, seems to confirm this observation. 

Let us finish by emphasizing that the nonlinear equation (28) is obtained 
in the closed form due to the large Nc limit assumption. In general, in the 
leading ln(l/x) approximation, a hierarchy of nonlinear equations is found 
by Balitsky [21]. These equations, written in a closed form by Weigert [32], 
are equivalent to the Color Glass Condensate renormalization group equa- 
tion [29], see [34] for details and more references. A similar equation to 
eq. (28) was also derived for diffractive DIS [35]. A different approach to 
unitarization at small x relies on the summation of interacting rcggcizcd 
gluons [36], forming a compound colour singlet state. The BFKL equation 
is obtained assuming two exchanged reggeized gluons. Recently, an energy 
spectrum of the scattering amplitude in high energy asymptotics with up to 
eight reggeized gluons was found in the large Nc limit [37] . For the discussion 
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of the relation between the approach based on multiple BFKL pomeron ex- 
changes, leading to eq. (28), and the approach based on exchanged reggeized 
gluons, see [38]. 

6. Conclusions and outlook 

We presented basic concepts of parton saturation. This phenomenon 
naturally appears in DIS at small x which is characterized by a strong rise 
of parton densities in the nucleon. However, the strong rise is tamed in 
accordance with unitarity of the description by the interactions between 
partons (mostly gluons). As a result, gluons form a coherent system char- 
acterized by the saturation scale, and the distribution of gluons with the 
transverse momentum below this scale no longer strongly rises. This effects 
arc important not only in ep DIS, but also in nucleus- nucleus collisions 
studied at RHIC [33]. 

Formulated in the dipole picture, the saturation effects impose strong 
bound on the form of the dipole-proton cross section for the dipole sizes big- 
ger than the inverse of the saturation scale. This feature was build in the 
phenomenological model [14] of the dipole cross section which successfully 
describes the inclusive and diffractive DIS data from HERA. In addition, 
based on the existence of the saturation scale, the new scaling ("geometric 
scaling") was proposed and successfully confronted with the data. The QCD 
based attempt to justify the main features of the phenomenological descrip- 
tion, which gives rise to a formulation of a nonlinear evolution equation, 
was presented. Based on the numerical solution to this equation we showed 
the main effect of saturation which suppresses diffusion into small values 
of transverse gluon momenta, present in the approach based on the linear 
BFKL equation. We also discussed the relevance of the impact parameter 
studies which touch the problem of proton confinement. Experimentally, 
this is done through the diffractive vector meson production. 

From the theoretical side, the future work will concentrate on detailed 
analysis of nonlinear equations describing parton saturation. Refraining 
from the large N(. limit, there are two equivalent formulations, the Colour 
Glass Condensate equation [29] and the equation derived by Weigert [32] 
which beautifully summarizes the infinite hierarchy of the equations ob- 
tained by Balitsky [21]. It would be good to have full handle over these 
equations. From the phenomenological side, more precise data awaited from 
HERA II will allow to test the parton saturation formulations in great detail. 
For this purpose, the phenomenological saturation model (14) was refined in 
order to include the DGLAP gluon evolution, expected at small r [39]. All 
important features of the saturation in the original model, however, were 
retained. 



18 



DIS3 PRINTED ON FEBRUARY 1, 2008 



The attractiveness of the presented description of saturation effects to 
large extent rehes on the dipole formulation (10), obtained in the leading 
ln(l/a;) approximation. How this picture changes when the next-to-leading 
order corrections to the formula (5) are computed, which introduce the qqg 
Fock state in the photon wave function, is not known yet. The computations 
are under way [40]. 
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